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, Abstract The paper deals with a third order semilinear equation which char- 

■ acterizes exponentially shaped Josephson junctions in superconductivity. The 

initial-boundary problem with Dirichlet conditions is analyzed. When the 
source term is a linear function, the problem is explicitly solved by means 
of a Fourier series with properties of rapid convergence. When F is nonlin- 
ear, appropriate estimates of this series allow to deduce a priori estimates, 
continuous dependence and asymptotic behaviour of the solution. 
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■ 1 Introduction 

a^: 

00 ■ We refer to the semilinear equation 

in 

fSj- (1) CeU = F{x,t,u) 

■ where is the third - order parabolic operator: 



X: 



(2) Ce = {d^x - \d..){£dt + 1) - dt{dt + a). 

The equation ([T]) characterizes the evolution of several dissipative models 



■ such as the motions of viscoelastic fluids or solids T-4j ; the sound propagation 
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in viscous gases [S]; the heat conduction at low temperature [5J[7] and the 
propagation of locaUzed magnetohydrodinamic models in plasma physics [5] . 
Moreover, it can also be referred to reaction diffusion systems [9 . 

As example of perturbed model of the phase evolution, we will consider the 
non linear phenomenon concerning the Josephson effects in superconductivity. 

More precisely, it ip = ip{x, t) is the phase difference in a rectangular junc- 
tion and 7 is the normalized current bias; when A = and F = sinip — 7, the 
equation ^ gives the well-known perturbed Sine-Gordon equation (PSGE) 



The terms £(fxxt and a (ft characterize the dissipative normal electron 
current flow respectively along and across the junction. They represent the 
perturbations with respect to the classic Sine Gordon equation [lOllllj . When 
the surface resistance is negligible, then e (< f ) is vanishing and a singular 
perturbation problem for the equation ([3]) could appear [12] . As for the coeffi- 
cient a of ([3]), it depends on the shunt conductance [13] and generally one has 
a < 1 [T^16j . However, if the resistance of the junction is so small as to short 
completely the capacitance, the case a > 1 arises [TUlfTT] . 

More recently, the case of the exponentially shaped Josephson junction 
(ESJJ) has been considered. The evolution of the phase inside this junction is 
described by the third order equation: 

(4) Sfxxt + ^xx - ^tt - e'^^Pxt - ^^x - aLpt = sin(^ - 7 

where A is a positive constant generally less than one [TSUin] while the terms 
A Lpx and A e ipxt represent the current due to the tapering. In particular X^px 
correspond to a geometrical force driving the fluxons from the wide edge to 
the narrow edge. [THlUn] . 

According to recent literature [l4l[T8l[20H23] , an exponentially shaped Joseph- 
son junction provides several advantages with respect to a rectangular junc- 
tion. For instance in [TB] it has been proved that in an ESJJ it is possible 
to obtain a voltage which is not chaotic anymore, but rather periodic exclud- 
ing, in this way, some among the possible causes of large spectral width. It is 
also proved that the problem of trapped flux can be avoided. Moreover, some 
devices as SQUIDs were built with exponentially tapered loop areas. [M] 

The analysis of many initial - boundary problems related to the PSGE 
(121) has been discussed in a lot of papers. In particular in [21], to deduce an 
exhaustive asymptotic analysis, the Green function of the linear operator 



(3) 



S^xxt 



+ Vxx - ^tt - OLipt = sin (/? - 7. 



(5) 



£ = dxx{edt + c^)-dt{dt + a) 



has been determined by Fourier series. By means of its properties an exponen- 
tial decrease of both linear and non linear solutions is deduced. 
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The aim of this paper is the analysis of the Dirichlet boundary value prob- 
lem related to the equation 

The Green function G of the linear strip problem is determined by Fourier 
series and properties of rapid convergence are established. So, when the source 
term is a linear function, then the explicit solution is obtained and an 
exponential decrease of the solution is deduced. 

When F is nonlinear, the problem is reduced to an integral equation with 
kernel G and an appropriate analysis implies results on the existence and 
uniqueness of the solution. Moreover, by means of suitable properties of G, 
a priori estimates, continuous dependence upon the data and asymptotic be- 
haviour of the solution are achieved, too. 

2 Statement of the problem and properties of the Green function 

Let /,T be arbitrary positive constants and let 

Ht = {{x, t) ■.0<x <l, 0<t<T}. 
The boundary value problem related to equation ([1]) is the following: 

{{dxx - A^o; ) {eut + u) - dt{ut +au) = F{x,t,u), {x,t) G fir, 
u{x,0) ^ ho{x), ut{x,0) ^ hi{x), xe[OJ], 
u(0,t)=0, w(/,i) = 0, 0<t<T. 

By Fourier method it is possible to determine the Green function of the 
linear operator C^- So, let 



(7) 7n = — , bn = {-fl + 5„ = -(a + e6„), = - &n 
and 

(8) Gn{t) = — e^f"* sinh{uJnt), 

by standard techniques, the Green function can be given the form: 

(9) G{x,t,^) ^ J e^'-" ^ Gnit) sinjn^ sin^/nX. 

This series is endowed of rapid convergence and it is exponentially vanishing 
as t tends to infinity. In fact, if we denote by 



(10) 
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and 



(11) PA 



the following theorem holds: 



q\ = , d = rmn{px,qx), 



Theorem 1 Whatever the constants a, e, A may be in the function G{x, ^, t) 
defined in (0i and all its time derivatives are continuous functions in Qt o,nd 
it results: 

(12) \G{x,U)\ < Me-'\ \^\<N,e-'', jSN 
where M,Nj are constants depending on a,X,e. 

Proof. Physical problems lead to consider ae < 1 and denoting by 



(13) 



^1,2 = 



I 



1/2 



let us assume that > 1- So, let fc be a positive constant less than one and 
let Ni^2,N^ be the lowest integers such that 



(14) 



4 (It Vl-ake) -e^kX^ 



1/2 



We start analysing the hyperbolic terms when n> N2- Letting 



(15) X„ = % < 1 



5n = ffn(-l + Vl - ^n), 



it is possible to prove that ipn < —. So it results: 

2.9ri 



(16) 



Furthermore, it is easily verified that for all n > ( > N2) it results < k 
and hence one has: 



(17) 



2f 



>5„ (l-fc)i/2 > (i_fc) 



1/2 
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Other terms can be treated similarly. For instance, as for circular terms, it 
can be proved that e~ ^" * < e""?^ *. 
In consequence estimate i holds V n > 1 . 

As for (HI]) 2, one has : 



b 2 

(18) 5„ - = < - + Vn > 1 

and by means of standard computations, (|T2t 9 can be deduced, too. 

It may be similarly proved that the theorem holds also when ae > 1 or 
when the conditions > 1 do not hold. 

Finally we notice that when are integers, the constants M and Nj in 
(IT^ could depend on t. 

As for the x-derivatives of Fourier series like ©, attention is needed to- 
wards convergence problems. For this, we will consider x-differentations of the 
operator {edt + 1)G instead of G and Gt- 

Theorem 2 Whatever a,e,X may be , the function G{x,^,t) defined in (0) 
is such that: 



(19) |a« {eGt + G)\ <A, e-'\ (i = 0,l,2) 

where 5 is defined in Ml]) and Ai {i = 0,1,2) are constants depending on 
a, e, A. 

Proof. As for the hyperbolic terms in G, it results: 



(20) eGt + G = e'^ V £_{[l-£(g„-c.„)]e""*- [l-e(5„ + c.„)]e-""*}. 
— ; t w„ 

n—l 

where according to (jl5|) . it results: 

I - e{gn - UJn) = 1 +£^n- 

So, by means of Taylor's formula, one has: 



(21) 1 - £ (g„ - a;„) = 1 - | 5„ X„ - ^g^ X^-^egn (1^- y)5/2 

Besides, it is possible to prove that Vrt > 1 one has: 



(22) X„ < ^, {c = l ^4 + AV e ^) 



2 
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and for all ji > c(l + ci)(ci > 0) it results: 



(23) 



(1 - -s""- 



ci(ci+2)]3/2 



So, taking into account that 



t: 9n^n — i ; 2 

2 aA + e 7„ 



ebl 



{ax +£7^)^' 
there exists a positive constant fci such that : 



(24) 



1 , aP ki . 



Estimates of theorem [T] together with (l24l) show that the series terms related 
to the operator eGt + G have order at least of n^"^. So it can be differentiated 
term by term with respect to x and the estimate (1191) can be deduced. 

As solution of the equation C/^v — we will mean a continuous function 
v{x,t) which has continuous the derivatives vt,Vtt, dx{svt + v), dxx{£Vt + v) 
and these derivatives verify the equation. 
So, we are able to prove the following theorem: 

Theorem 3 The function G{x, t) defined in (0) is a solution of the equation 



(25) 



CeG = (dxx - Xdx)ieGt + G) - dtiGt + aG) = 0. 



Proof. The uniform convergence proved in theorems [T][2] allows to deduce 



that: 
(26) 



dG 2e^ °° 



dt I 



[6„(eff„-l)G„-e6„e cosh w„t] sin 7„^ sin 7„a 



2 °° r 

(27) dx{edt + \)G^-e^Y.\ 

n=l ^ 



[l-e gn ) Gn+ 



Xe 



_|_ _g 9r,t coshWni 



sin7„^sin7„a;+ 
+ [Gn In (1 - sgn) + e7„e"^"* cosh w„i] sin7„,^cos7„x 
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Moreover, beeing 



°° f \2 
(28) drAedt + l)G = ye^ [(^ 5n - l) (&n - y) G„- 

n— 1 ^ 



+ £7^)e"^"* coshwni] sin 7„^ sin 7„x- 



+ [G„A7„(1 - £g„) + £A7„e cosha;„t] sin7„^cos7„a; 
3 can be deduced. 



3 Properties of the convolution 

To achieve the sohition of the strip problem ([6]) , the convolution of the function 
G with the data must be analysed. For this, let h{x) be a continuous function 
on (0, 1) and let: 

(29) u^x^t) = [ h{OG{x,^,t) 







(30) ul^ix.t) ^ {dt+a + e\dx- edxx)uh{x,t). 

The following theorems hold: 

Theorem 4 If the data h{x) is a C^{0,1) function, then Uh defined by \2S\) is 
a solution of the equation = Q and it results: 



(31) \\muh{x,t)=Q lim dtUhix, t) — h(x), 

uniformly for all x G [0, /] . 

Proof. The absolute convergence of with its partial derivatives is proved 
by means of theorems [T] and [5] and continuity of function h{x). So, since ([25]) 
CgUh = is verified, while theorem [1] and hypotheses on h{x) imply (j3ip i . 

More beeing: 



dG 2 d ^ dGn cos7„^ . 
ot n ot ot n 



and 
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(33) 



dtUh 



2^ dGn . w sin7„a; 



TT — at 

n—1 



2 ^ dGn , C0S7n^ 



E 



./n ^ dt 

n—1 



■MO' 



denoting by ri{x) the Heaviside function, it results: 



(34) \imdtUh = j[hil)-h{0)]+h{0)-[ h'iO[vii-x) + j-l]d^ = h{x). 



Theorem 5 Let h{x) be a C3(0,/) function such that /iW(0) = h'^'^l) = 
(i = 1,2,3). Then defined in h3U\l is a solution of the equation £e = 
and it results: 



lim u^(x, t) = ft. lim 9tiiJ^(x, t) = 



(35) 

uniformly for all x G (0, 1). 

Proof. Properties of h{x) assure that: 



(36) 



{\dx ~ dxx)uh {x,t) 



■jiXdx-d^x) e^ VG„(t) / h'iO 
' , 1 Jo 



^ — sin7„x 



e 2 



E 

ra=l 



G„(i) hi 



h (C) 



sin7„^ 
7^ 



= - U,^n{x,t) + — Uhix,t) 

So, since theorem [31 £eujj = is verified. Moreover, beeing: 



(37) dtul = {dxx - >^dx)uh 

([55)1 imphes ([55)) 9. too. Finally, owing to (PT|) and ([55t , one obtains: 



(38) 



lim ul — lim 



h{x) 
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4 Solution of the linear problem 

Let us consider the homogeneous case. From theorems 21 [S] the following result 
is obtained: 

Theorem 6 When F ~ and the initial data hi{x), and ho^x) verify the 
hypotheses of theorems \S[ then the function: 

(39) u{x, t) = Uh^ + {dt+a + eXdx - ed^x) Uho 
represents a solution of the homogeneous strip problem (0). 

Otherwise, when F = f{x,t) , let consider 

(40) Uf{x,t) = - / dT [ /(e,r) G{x,^,t-r) d^. 

JQ Jo 

Standard computations lead to consider at first the problem ([B]) with go — 
gi = 0. For this the following theorem is proved: 

Theorem 7 If the function f{x,t) is a continuous function in Qt with con- 
tinuous derivative with respect to x, then the function u / represents a solution 
of the nonhomogeneous strip problem. 

Proof. Since (pijl i it results: 

(41) dtuf{x,t) = f dT f f{£,,T) Gt{x,^,t-T) d^ 

Jo JQ 

and as proved in theorem 21 one obtains: 

(42) \imdtUfix,t) = f{x,t). 
Hence, one has: 

(43) =f{x,t)+ f dT f f{^,T) Gu{x,^,t-T) d^ 

Jo Jo 

and theorem [31 assures that C^Uf = f{x,t) . 

Furthermore, owing to (j40D - (PT|) and estimates (fT^ . if Bi (i=l,2) are two 
positive constants, it results: 

(44) <Bi(l-e-**); < ^2(1 - e"^*) 
from which initial homogeneous conditions follow. 

The uniqueness is a consequence of the energy-method and we have: 
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Theorem 8 When the source term f{x,t) satisfies theorem^ and the initial 
data (/io,/ii) satisfy theorem\^ then the function 

(45) u{x, t) = u/ij + (9t + a + e Xd^ - ed^x) w/io + 

is the unique solution of the linear non-homogeneous strip problem 

5 Solution of the non hnear problem 

Let us consider now the non linear problem: 

{{da:x - A9j; ) [eut + u) - dt{ut + au) ^ F{x,t, u), {x,t) e Ht, 
u{x, 0) = ho{x), utix, 0) = hi{x), x G [0, 1], 
w(0,t) = 0, u{l,t)^0, 0<t<T, 

As for the data F and hi (x) {i = 0, 1) we shall admit: 

Assumption 5.1 The functions hi{x){i — 0,1) are continuously differ- 
entiable and bounded together with h'-^{x) and h\^^ {k — 1,2). The function 
F{ x,t,u) is defined and continuous on the set 

(47) Dt = { {x, t, u) : (x, t) e Qt , -oo < u < oo} 

and more it is uniformly Lipschitz continuous in {x, t, u) for each compact 
subset of f^T ■ Besides, F is bounded for bounded u and there exists a constant 
Cp such that the estimate 



(48) \F{x,t,ui) - F{x,t,U2)\ < Cf |mi - ^2 | 

holds for all (ui, U2). 

When the problem admits a solution u then, properties of G and the 
assumptions 5.1, assure that u must satisfy the integral equation 

(49) 

u{x,t)= [ hi{OG{x,i,t)di + {dt + a + e\dx~edxx) f ho{^)G{x,^,t)d^+ 

JQ Jo 

+ f dT f G{x,^,t-T)F{^,T,u{^,T))d^, 

Jo Jo 
and it is possible to prove that: 

Theorem 9 The non linear problem |^6p admits a unique solution if and only 
if the integral equation J^ffj ) has a unique solution which is continuous on Qt- 
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Moreover, let ||w||t = sup i) | and let Bt denote the Banach space 

(50) Bt = {vix,t) : w e C( J7t), ||w||t < oo }. 

By means of standard methods related to integral equations it is possible to 
prove that the mapping ip defined by (|49l) is a contraction of Bt in Bt and so 
it admits a unique fixed point u{x,t). In consequence the following theorem 
holds: 

Theorem 10 When the initial data hi [i = 0, 1) and the source term F verify 
the assumption 5.1, then the problem admits a unique regular solution. 

6 Applications 

All these results allow us to obtain continuous dependence upon the data, a 
priori estimates of the solution and asymptotic properties. 
According to assumption 5.1, let 



= sup \h,{x)\ , (« = 0, 1) 


K 


— sup 




(OJ) 




(0,0 





\\u\\rp = sup f)| , = sup\F{x,t,u)\ . 

So, by means of the following theorem the dependence upon the data can 
be proved: 

Theorem 11 Let ui, U2 be two solutions of the problem related to the data 
{ho, hi, Fi) and (70, 71 F2) which satisfy the assumption 5.1. Then, there 
exists a positive constant C such that 

\\ui — U2 \\t < C sup \ho — 7o| + C sup — 7i| + C sup \Fi{x, t, u) — F2{x, t,u)\ , 

(2t (2t Dt 

where C depends on Cp, T and on the parameters a, £, A. 

The integral equation and the properties proved for Green Function G 
imply a priori estimates, too. 

Theorem 12 When the data {ho,hi, F) of the problem |^6p verify the as- 
sumption 5.1, then the following estimate holds: 

(51) \\u{x,t)\\^ < i (1-e-**) ||F|| + K [ ||;.i|| + ll/ioll + \\ho\\ le"^' 
where the constants S- defined in ill}) - and K depend on a,e, A. 
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As for the asymptotic properties, obviously the behaviour of the solution de- 
pends upon the shape of the source term. 
For instance, in the linear case one has: 

Theorem 13 When the source term f{x,t) satisfies the condition: 

(52) < Ce^™* {C, m = const >0), 
one has; 

(53) |u(a;,t)| < fce~™* m*^min{S,m}, k = const. 

An exponentially decreasing behaviour is also possible in the non linear 
case. In fact, according to [29j . let us consider a normed space where 

(54) ||M(x,t)|| = max |M(x,i)| 

xe{o,i) 

is such that 

(55) \\u{x,t)\\ < P e-'' 

beeing [3 a positive constant and 6 is defined in ([TT]) . Furthermore, let us 
introduce the following definition ,25] : 

Definition 7.1 When the function F is such that \F{x, t,u)\ < 7 ||u|| e""^*, 
then F is an exponential Lipschitz function. 

So the following theorem can be proved: 

Theorem 14 // the non linear source F is an exponential Lipschitz function, 
then the solution of the semilinear problem 1^6] ) vanishes as follows: 

(56) \u{x,t)\ < Ki e"** 
where Ki is a positive constant depending on a, e, A. 

Since | sin u\ < a similar behaviour is also verified for the model of super- 
conductivity when F(x, t, u) = sin u. 

Acknowledgements This paper has been performed under the auspices of 
G.N.F.M. of I.N.D.A.M. 
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